Critical phenomenon of the single crystal helimagnet MnSi 
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Abstract 


The critical behavior of the single crystal helimagnet MnSi is investigated by means of bulk 


dc-magnetization at the boundary between the conical state and paramagnetic phase. We obtain 


the critical exponents (3 = 0.242 + 0.006, y = 0.915 + 0.003, and 6 = 4.734 + 0.006), where the 


self-consistency and reliability are verified by the Widom scaling law and the scaling equation. 


The critical exponents of MnSi belong to the universality class of tricritical mean-field theory, 


which unambiguously indicates a tricritical phenomenon at the boundary between the first-order 


phase transition and the second-order one induced by the external magnetic field. The tricritical 


point (TCP) is determined as Hrcp = 3200 Oe at the critical temperature, consistent with the 


previous report [Phys. Rev. Lett. 110, 177207 (2013)]. The critical behavior suggests a long-range 


magnetic coupling with the exchange distance decaying as J(r) ~ r743 in MnSi. 
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I. INTRODUCTION 


For decades, the B20 transition metal compound MnSi has always been fascinating due 
to its weak itinerant ferromagnetism [1], non-Fermi-liquid behavior |2} [3], and magnetic 
quantum phase transition |4|. Recently, there is a renewed interest in MnSi because of 
the discovery of the skyrmion phase |5} [6] 17, 8) 9) (11), which has great potential for 
new spin-based applications on account of the topological property |5, [6], nanometric size 
13], and current-driven motion |[14, 117]. The crystal structure of MnSi is a non- 
centrosymmetric cubic cell with space group P213, which results in a weak Dzyaloshinskii- 
Moriya (DM) spin-orbit coupling due to the lack of inversion symmetry. The DM spin-orbit 
coupling competes with the much stronger ferromagnetic exchange, which finally causes a 
long modulation period of a helimagnetic ground state |5, 6,118]. When an external magnetic 
field is applied above a threshold value, the helimagnetic state transforms into a conical phase 
[19]. The chiral skyrmion phase emerges just below Tọ in a narrow magnetic field range 
due to a unique stabilization provided by the DM interactions [5]. The ordering processes 
in cubic helimagnets are characterized by the formation of specific precursor states |18} 20]. 
It is demonstrated that a disordered Skyrmion phase appears in a very narrow temperature 
range of ~ 1 K above the helical phase in MnSi [23]. 

Recently, the critical phenomena of some helimagnets with a skyrmion phase have been 
analyzed intensively to investigate the detailed ferromagnetic exchange. For example, it has 
been demonstrated that the near-neighbor isotropic exchange of short-range type is domi- 
nant in Feo gCoo2Si [24]. In addition, another multiferroic insulating helimagnet Cu2OSeOs, 
which also exhibits a skyrmion state, is consistent with the classical three-dimensional (3D) 
Heisenberg model with short-range interaction |25|. However, for MnSi, the magnetic ex- 
change is expected to be different from that in Fe;_,Co,Si or CugOSeO3 because a long-range 
coupling as described by the mean-field theory has been predicted by the theory of Brazovskii 
[26]. According to the Brazovskii scenario, a first-order phase transition from paramagnetic 
to helimagnetic ordering was induced by the critical fluctuation under zero-field |26|, which is 
supported experimentally by the neutron scattering studies |27|, the renormalization group 
theory and Monte Carlo simulations |28|. The weak first-order phase transition can 
be suppressed into a second-order one by the external magnetic field, prospectively yielding 


tricritical behavior [29,30]. However, the tricritical phenomenon of this system has not been 


studied in detail. In this work, the critical behavior of the single crystal MnSi is investi- 
gated around To by means of bulk dc-magnetization. We find that the critical exponents 
of MnSi satisfy the universality class of the tricritical mean-field theory, which unambigu- 
ously indicates a tricritical phenomenon at the boundary of the first-order transition and the 
second-order one. The tricritical point is determined to be Hrcp œ% 3200 Oe at To, which 
agrees with the results of specific heat experiments [80]. Moreover, the critical behavior 
indicates a long-range magnetic coupling with the magnetic exchange distance decaying as 


Jr) = r74. 


II. EXPERIMENT 


A single crystal sample of MnSi was prepared by the Czochralski method [31]. The 
structure and phase purity were checked with a Rigaku-TTR3 X-ray diffractometer. The 
chemical compositions were carefully determined by energy dispersive x-ray spectrometry 
(EDXS), which shows the proportion of Mn: Si ~ 49: 51. The magnetization was measured 
using a Quantum Design vibrating sample magnetometer (SQUID-VSM). The no-overshoot 
mode was applied to ensure a precise magnetic field. In addition, the field was relaxed for 
2 minutes before the data collection. For the measurement the sample was processed into 
slender ellipsoid shape, and the magnetic field was applied along the longest axis to minimize 
the demagnetizating field. To make sure each curve was initially magnetized, the isothermal 
magnetization was performed after the sample was heated well above Teo for 10 minutes 
and then cooled under zero field to the target temperatures. The magnetic background was 
carefully subtracted. The applied magnetic field H, has been corrected into the internal field 
as H = Ha — NM (where M is the measured magnetization and N is the demagnetization 


factor) [32]. The corrected H was used for the analysis of critical behavior. 


II. RESULTS AND DISCUSSION 


Figure [I| (a) shows the temperature dependence of magnetization M(T) for MnSi un- 
der zero-field-cooling (ZFC) and field-cooling (FC) with H = 100 Oe (left coordinate). 
A paramagnetic-helimagnetic (PM-HM) transition occurs at To = 30.5 K determined by 
dM/dT (right coordinate), which agrees with the expected mean-field transition tempera- 


ture Tyr % 30.5 K 133]. A narrow dip appears at ~ 29.7 K just below To (see the inset 
of Fig. | (a)), which is a clue to the first-order phase transition. Figure [I] (b) gives the 
isothermal magnetization M(H) at 4 K, which exhibits a magnetic ordering behavior with 
the saturation field Hs ~ 6000 Oe. The inset of Fig. [1] (b) plots the magnification of M(H) 
in the lower field region, which indicates almost no coercive force for MnSi. 

As is well known, the critical behavior for a second-order phase transition can be inves- 
tigated in detail through a series of critical exponents. In the vicinity of a second-order 
phase transition, the divergence of correlation length € = &|(T — Tc¢)/To|~” leads to uni- 
versal scaling laws for the spontaneous magnetization Ms and initial susceptibility xo. The 


mathematical definitions of the exponents from magnetization are described as [84 35]: 


Ms(T) = Mo(-e)’,e < 0,T < Te (1) 
Xo (T) = (ho/Mo)e",e > 0,T > To (2) 
M = DH!’ e=0,T=Tec (3) 


where € = (T — Tc)/Tc is the reduced temperature; Mo/ho and D are critical amplitudes. 
The parameters 8 (associated with Ms), y (associated with xo), and ô (associated with To) 
are the critical exponents. Generally, the critical exponents should follow the Arrott-Noakes 


equation of state |36): 
(H/M)") = (T — Te) /To + (M/Mı)"? (4) 


Therefore, the critical exponents 3 and y can be obtained by fitting the Ms(T) and xg (T) 
curves using the modified Arrott plot of MV/P vs (H/M)'/7. Meanwhile, the third exponent 
ô can be generated directly by Eq. (3). 

For the analysis of the critical behavior, we measure the initial isothermal M(H) curves 
around Tc, which are shown in Fig. [2] (a). Previous theoretical investigation suggested that 
magnetic coupling in MnSi should be described using the mean-field theory |26). In the 
mean-field theory where 8 = 0.5 and y = 1.0 [87], we generate the Arrott plot of M? vs 
H/M, as shown in Fig. B| (b). Generally, M? vs H/M should be a series of parallel straight 


lines in the higher field range. The lower field data, which mainly represent the arrangement 


of magnetic domains, should be excluded in the fitting process [38]. The line of M? vs H/M 
at Tc should pass through the origin [39]. According to the criterion proposed by Banerjee, 
the order of the magnetic transition can be determined from the slope of the straight line: 
a positive slope corresponds to the second-order transition while the negative to the first- 
order one |40]. Apparently, the positive slopes of the present M? vs H/M curves imply that 
the magnetic transition of MnSi in the higher field region is a second-order one. However, 
all the curves in the Arrott plot are nonlinear and show downward curvature even in the 
high field region, in agreement with previous report [41]. This suggests that the mean-field 
theory with G= 0.5 and y= 1.0 is not appropriate for MnSi according to the Arrott-Noakes 
equation of state [36]. Hence, a modified Arrott plot should be employed to obtain the 
critical exponents. 

Three kinds of possible exponents belonging to the 3D-Heisenberg model (G= 0.365, 
y=1.336), 3D-Ising model (G= 0.325, y=1.24), and tricritical mean-field model (G= 0.25, 
y=1.0) are used to make a modified Arrott plot, as shown in Fig. B| (a), (b), and (c). All 
three models exhibit quasi-straight lines in the high field region [42,143] |44|. To determine an 
appropriate model, the modified Arrott plots should be a series of parallel lines in the high 
field region with the same slope (where the slope is defined as S(T) = dM'/8/d(H/M)'/7). 
The normalized slope (VS) is defined as NS = S(T)/S(Tc), which enables one to identify 
the most suitable model by comparing NS with the ideal value of ’1’ |42|. Plots of NS vs 
T for the three different models are shown in Fig. [B] (d), which indicates that the tricritical 
mean-field model is the best one to describe the critical behavior of MnSi. 

The linear extrapolation from the high field region to the intercepts with the axes M1/? 
and (H/M)'/7 yields reliable values of spontaneous magnetization Ms(T,0) and inverse 
initial susceptibility xg (T,0), which are plotted as a function of temperature in Fig. 4 (a). 
By fitting to Eq. and Eq. (2), one obtains 3 = 0.242+0.006 with To = 30.38 + 0.04 and 
y = 0.915 + 0.003 with To = 30.69 + 0.07. The critical temperature To from the modified 
Arrott plot is very close to that obtained from the M(T) curve. Figure [4| (b) shows the 
isothermal magnetization M(H) at To = 30.5 K, where the inset plots are on a lg — lg scale. 
We determine that 6 = 4.734 + 0.006 in the high field region (H > Hs). It can be seen 
that the M(H) at To exhibits a straight line on a lg — lg scale for H > Hg. According to 


statistical theory, these critical exponents should fulfill the Widom scaling relation [45]: 


et 
er, (5) 


As a result, 6 = 4.781 + 0.007 is calculated following the Widom scaling relation based 
on the obtained 3 and y, which agrees well with the results from the experimental critical 
isothermal analysis. The self-consistency of these critical exponents demonstrates that they 
are reliable and unambiguous. In order to further testify the convergence of the critical 


exponents, the effective exponents bepp and Yerr for MnSi are obtained as [46]: 


= 

Berro = TEP urn) = Se (6) 
The effective exponents d.rp and Yerp as a function of reduced temperature £ for MnSi are 
plotted in Fig. It can be seen that Ger, and yer are convergent when the temperature 
approaches To. 

The obtained critical exponents of MnSi, as well as those from different theoretical models 
and related materials, are listed in Table I for comparison. It can be seen that the expo- 
nents of MnSi are very close to those predicted by the tricritical mean-field theory, which 
unambiguously indicates a tricritical phenomenon. The tricritical phenomenon usually oc- 
curs at the boundary between a first-order phase transition and a second-order one, where 
the three lines of critical transitions meet at a single point: the so called tricritical point 
(TCP) [47 [49|. In the manganese oxides, the tricritical phenomenon usually occurs 
when the first-order phase transition disappears and the second-order transition is induced 
by proper doping effects |47, |48) [49]. In the present case of MnSi, a similar phase transition 
from first-order to second-order also occurs. It should be noted that an isotropic critical 
fluctuation induced first-order phase transition happens to MnSi under zero field, which can 
be easily suppressed by an external magnetic field 28]. Thus, in this case, the 
tricritical phenomenon is expected at the boundary of these two phases. The critical point 
appears when the first-order phase transition is suppressed into a second-order one by an 
external field. The existence of a field-induced tricritical phenomenon was also pointed out 
by A. Bauer et al. in terms of the high-precision measurement of the specific heat [80]. In 
fact, high pressures can also induced a tricritical behavior in MnSi with P;, & 3.5 GPa and 
Tr ~ 25 K, where the pressure modulates the orders of the phase transitions 52]. 


The critical exponents obtained here in the single crystal sample are slightly different from 


those in the polycrystalline one, where the exponent 6 was calculated using the Widom 
scaling relation rather than obtained from the experiment [53]. As can be seen in Table I, 
the exponents obtained in the polycrystalline sample cannot be categorized into any uni- 
versality classes, which may be attributed to the magnetic scattering from lots of crystal 
boundaries and defects. As we know, when the temperature approaches to the critical point, 
the correlation length € exceeds the grain size, where the critical exponents are affected by 
the crystal boundaries and defects |54|. Thus, the intrinsic magnetic coupling is masked in 
an inhomogeneous polycrystalline sample, especially for the long-range interaction [54]. 
The agreement of MnSi with the tricritical mean-field model indicates a long-range mag- 
netic coupling which is revealed by the exponent y |55|. As we know, for a homogeneous 
magnet, the universality class of the magnetic phase transition depends on the exchange 
distance J(r). M. E. Fisher et al. theoretically treated this kind of magnetic ordering as an 
attractive interaction of spins, where a renormalization group theory analysis suggests the 


interaction decays with distance r as |55) 56]: 
Kr) art) (7) 


were d is the spatial dimensionality and ø is a positive constant. Moreover, there is [32) 55} 
57]: 

in g so a nr 8) 2 
where Ao = (ø — $) and G(Ż) = )?, 


2 


2G(4)(7n + 20) 
(n — 4)(n +8) 


(8) 


n is the spin dimensionality. For a three 
dimensional material (d = 3), we have J(r -8+0), When o > 2, the Heisenberg model 
(8 = 0.365, y = 1.386 and 6 = 4.8) is valid for the three dimension isotropic magnet, where 
J(r) decreases faster than r~°. When o < 3/2, conditions for the the mean-field model 
(8 = 0.5, y = 1.0 and 6 = 3.0) are satisfied, expecting that J(r) decreases slower than 
r745, From Eq. |8) we see that o = 1.329 + 0.008 for MnSi, thus confirming the long-range 
magnetic coupling for o < 3/2. Subsequently, we find that the exchange distance decays 
as J(r) œ~ r743, which indicates slow decay of the magnetic coupling. The decay distant of 
MnSi suggests that the spatial magnetic interaction is relatively extensive. Moreover, we 
get the correlation length critical exponent v = 0.688 + 0.009 (where v = y/o), which is 
close to 0.62 + 1 obtained by small-angle polarized neutron scattering (SAPNS) [58]. 
Finally, these critical exponents should follow the scaling equation. According to the 


prediction of scaling equation, in the asymptotic critical region, the magnetic equation can 


7 


be written as [34]: 
M(H, «) = ef fa (H/e°*7) (9) 


where f+} are regular functions denoted as f} for T > To and f_ for T < Tec. Defining 


the renormalized magnetization as m = e”?M(H,e), and the renormalized field as h = 


He-(+N), the scaling equation indicates that m vs h forms two universal curves for T > To 


and T < To, respectively |49) 54]. Based on the scaling equation [m = f+(h)], the isothermal 


magnetization around the critical temperatures for MnSi is plotted in Fig. [6] (a), where all 
experimental data in the higher field region collapse onto two universal curves. In addition, 
the experimental data in low field region above T¢ collapse well. However, the lower-field 
experimental data below Tç cannot collapse onto one curve, as shown in the inset of Fig. 
6l (a). The m? vs h/m curves [Fig. [6] (b)] apparently exhibit a vertical line in the lower 
field region. As we know, MnSi exhibits a first-order phase transition under a low external 
magnetic field, which is suppressed into a second-order one with a strong external field. 
The TCP of magnetic field (Hrcp) lies at the boundary between the first-order and the 
second-order transition. Above Hrgp in the second-order transition region, the M — T — H 
curves should collapse onto a universal curve. However, below Hrcp in the first-order 
transition region, the scaling equation is invalid. Thus, the turning point of the external 
field He on the m? vs h/m curve reflects the boundary between the first-order transition 
and the second-order one, as shown in Fig. [6| (b). The Hc vs T is plotted in the inset of 
[6] (b), which decreases almost linearly with increasing temperature when approaching To. 
The point where the Hg and To meet is the Hrcp. The linear fitting of the Hco(T) yields 
Hrcp ~ 3200 Oe at To, which is very close to 3400 Oe obtained by specific heat experiments 
30). 

As mentioned above, the critical exponents of MnSi indicate a long-range magnetic 
coupling. However, as can be seen in Table I, the critical exponents of Feo.gCop.25i and 
Cuz,0SeO3, which also exhibit a helimagnetic and skyrmion phase transition with similar 
crystal symmetry, are close to the universality class of the 3D-Heisenberg model [24, 25]. 
That is to say, the micro-magnetic interactions in Feg..Cog.2Si and CugOSeQO3 are of the 
short-range type. The different critical behaviors of MnSi from Feg.sC00.2Si or CugOSeO3 
indicate that the micro-magnetic interactions are intrinsically different. The micro-magnetic 
interaction may correlate with the transport behaviors of the system. The weak itinerant 


ferromagnetism in MnSi may probably be related to the long-range magnetic coupling since 


it has been proved that the transport behavior is connected to the miro-magnetic interaction 


in these systems [24]. 


IV. CONCLUSION 


In summary, the critical behavior of the single crystal helimagnet MnSi has been investi- 
gated around Tg. We obtained the critical exponents 8 = 0.242 + 0.006, y = 0.915 + 0.003, 
and 6 = 4.734 + 0.006, which are verified with the Widom scaling law and the scaling equa- 
tion. The critical exponents of MnSi belong to the universality class of tricritical mean-field 
theory, which unambiguously indicates a tricritical phenomenon. It is suggested that the 
field-controlled phase transition from the critical fluctuation induced first-order to second- 
order should be responsible for the appearance of the tricritical phenomenon in MnSi. Based 
on these critical exponents, the critical exponent for the correlation length is deduced to 
be v = 0.688 + 0.009 (€ = &|(T — Te)/Te|””). The tricritical point is determined as 
Hrop = 3200 Oe at Ty. The exchange distance is obtained to decay as J(r) ~ r743, which 


suggests a long-range magnetic coupling in MnSi. 


V. ACKNOWLEDGEMENTS 


This work was supported by the State Key Project of Fundamental Research of China 
through Grant Nos. 2010CB923403 and 2011CBA00111, the National Natural Science Foun- 
dation of China (Grant Nos. U1332140, U1232142, 11004196, 11474290, 11104281, and 
11204288). 


1] N. Manyala, Y. Sidis, J. F. DiTusa, G. Aeppli, D. P. Young, Z. Fisk, Nature (London) 404, 
581 (2000). 

2] C. Pfleiderer, S. R. Julian, G. G. Lonzarich, Nature (London) 414, 427 (2001). 

H. Watanabe, S. A. Parameswaran, S. Raghu, A. Vishwanath, Phys. Rev. B 90, 045145 (2014). 


A WwW 


C. Pfleiderer, G. J. McMullan, S. R. Julian, G. G. Lonzarich, Phys. Rev. B 55, 8330 (1997). 


5] U. K. Roßler, A. N. Bogdanov, C. Pfleiderer, Nature (London) 442, 797 (2006). 


[6] 


13 
14 


15 


16 


17 


18 
19 
20 


21 


22 


S. Muuhlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, P. Boni, 
Science 323, 915 (2009). 

W. Munzer, A. Neubauer, T. Adams, S. Muhlbauer, C. Franz, F. Jonietz, R. Georgii, P. Boni, 
B. Pedersen, M. Schmidt, A. Rosch, C. Pfleiderer, Phys. Rev. B 81, 041203 (2010). 

X. Yu, Y. Onose, N. Kanazawa, J. Park, J. Han, Y. Matsui, N. Nagaosa, Y. Tokura, Nature 
(London) 465, 901 (2010). 

S. Seki, X. Z. Yu, S. Ishiwata, Y. Tokura, Science, 336, 198 (2012). 

H. F. Du, W. Ning, M. L. Tian, Y. H. Zhang, Phys. Rev. B 87, 014401 (2013). 

A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. G. Niklowitz, P. Boni, Phys. Rev. 
Lett. 102, 186602 (2009). 

H. F. Du, J. P. DeGrave, F. Xue, D. Liang, W. Ning, J. Y. Yang, M. L. Tian, Y. H. Zhang, 
S. Jin, Nano Lett. 14, 2026 (2014). 

N. Nagaosa and Y. Tokura, Nat. Nanotechnol. 8, 899 (2013). 

F. Jonietz, S. Muhlbauer, C. Pfleiderer, A. Neubauer, W. Munzer, A. Bauer, T. Adams, R. 
Georgii, P. Boni, R. A. Duine, K. Everschor, M. Garst, A. Rosch, Science 330, 1648 (2010). 
X. Z. Yu, N. Kanazawa, W. Z. Zhang, T. Nagai, T. Hara, K. Kimoto, Y. Matsui, Y. Onose, 
Y. Tokura, Nat. Commun. 3, 988 (2012). 

A. Fert, V. Cros, J. Sampaio, Nat. Nanotechnol. 8, 152 (2013). 

J. S. White, K. Prsa, P. Huang, A. A. Omrani, I. Zivkovic, M. Bartkowiak, H. Berger, A. 
Magrez, J. L. Gavilano, G. Nagy, J. Zang, H. M. Ronnow, Phys. Rev. Lett. 113, 107203 
(2014). 

P. Bak and M. H. Jensen, J. Phys. C: Solid St. Phys. 13, L881 (1980). 

A. Bauer and C. Pfleiderer, Phys. Rev. B 85, 214418 (2012). 

H. Wilhelm, M. Baenitz, M. Schmidt, C. Naylor, R. Lortz, U. K. Ro@ler, A. A. Leonov, A. 
N. Bogdanov, J. Phys.: Condens. Matter 24, 294204 (2012). 

C. Pappas, E. Lelievre-Berna, P. Falus, P. M. Bentley, E. Moskvin, S. Grigoriev, P. Fouquet, 
B. Farago, Phys. Rev. Lett. 102, (2009) 197202. 

S. V. Grigoriev, N. M. Potapova, S. A. Siegfried, V. A. Dyadkin, E.V. Moskvin, V. 
Dmitriev, D. Menzel, C. D. Dewhurst, D. Chernyshov, R. A. Sadykov, L. N. Fomicheva, 
A.V. Tsvyashchenko, Phys. Rev. Lett. 110, 207201 (2013). 


[23] S. V. Grigoriev, S. V. Maleyev, E. V. Moskvin, V. A. Dyadkin, P. Fouquet, H. Eckerlebe, 


10 


[24] 
[25] 


26 
27 


28 
29 
30 
31 


32 
33 


34 


35 
36 
37 
38 
39 
40 
41 
42 


43 


44 


[45] 
[46] 


Phys. Rev. B 81, 144413 (2010). 

W. J. Jiang, X. Z. Zhou, G. Williams, Phys. Rev. B 82, 144424 (2010). 

I. Zivkovic, J. S. White, H. M. Ronnow, K. Prsa, H. Berger, Phys. Rev. B 89, 060401(R) 
(2014). 

S. A. Brazovskii, Sov. Phys. JETP 41, 85 (1975). 

M. Janoschek, M. Garst, A. Bauer, P. Krautscheid, R. Georgii, P. Boni, C. Pfleiderer, Phys. 
Rev. B 87, 134407 (2013). 

S. Buhrandt and L. Fritz, Phys. Rev. B 88, 195137 (2013). 

K. Huang, Statistical Mechanics, 2nd ed. (Wiley, New York, 1987). 

A. Bauer, M. Garst, C. Pfleiderer, Phys. Rev. Lett. 110, 177207 (2013). 

V. A. Dyadkin, S. V. Grigoriev, D. Menzel, D. Chernyshov, V. Dmitriev, J. Schoenes, S. V. 
Maleyev, E. V. Moskvin, H. Eckerlebe, Phys. Rev. B 84, 014435 (2011). 

A. K. Pramanik and A. Banerjee, Phys. Rev. B 79, 214426 (2009). 

J. Kindervater, W. Hauler, M. Janoschek, C. Pfleiderer, P. Boni, M. Garst, Phys. Rev. B 
89, 180408(R) (2014). 

H. E. Stanley, Introduction to Phase Transitions and Critical Phenomena (Oxford University 
Press, London, 1971). 

M. E. Fisher, Rep. Prog. Phys. 30, 615 (1967). 

A. Arrott and J. Noakes, Phys. Rev. Lett. 19, 786 (1967). 

S. N. Kaul, J. Magn. Magn. Mater. 53, 5 (1985). 

N. Khan, P. Mandal, K. Mydeen, D. Prabhakaran, Phys. Rev. B 85, 214419 (2012). 

A. Arrott, Phys. Rev. 108, 1394 (1957). 

S. K. Banerjee, Phys. Lett. 12, 16 (1964). 

D. Bloch, J. Voiron, V. Jaccarino, J. H. Wernick, Phys. Lett. 51A, 259 (1975). 

J. Y. Fan, L. S. Ling, B. Hong, L. Zhang, L. Pi, Y. H. Zhang, Phys. Rev. B 81, 144426 (2010). 
L. Zhang, B. S. Wang, Y. P. Sun, P. Tong, J. Y. Fan, C. J. Zhang, L. Pi, Y. H. Zhang, Phys. 
Rev. B 85, 104419 (2012). 

L. Zhang, J. Y. Fan, L. Li, R. W. Li, L. S. Ling, Z. Qu, W. Tong, S. Tan, Y. H. Zhang, 
Europhys. Lett. 91, 57001 (2010). 

L. P. Kadanoff, Physics 2, 263 (1966). 

A. Perumal, V. Srinivas, V. V. Rao, R. A. Dunlap, Phys. Rev. Lett. 91, 137202 (2003). 


11 


[47] 


48 
49 


50 
51 
52 


[53] 
[54] 


55 
56 


57 
58 


D. Kim, B. Revaz, B. L. Zink, F. Hellman, J. J. Rhyne, J. F. Mitchell, Phys. Rev. Lett. 89, 
227202 (2002). 

J. Yang, Y. P. Lee, Y. Li, Phys. Rev. B 76, 054442 (2007). 

M. H. Phan, V. Franco, N. Bingham, H. Srikanth, N. Hur, S. Yu, J. Alloys Compd. 508, 238 
(2010). 

A. E. Petrova, V. N. Krasnorussky, J. Sarrao, S. M. Stishov, Phys. Rev. B 73, 052409 (2006). 
A. E. Petrova, V. N. Krasnorussky, J. Sarrao, S. M. Stishov, JETP 102, 636 (2006). 

M. Otero-Leal, F. Rivadulla, S. S. Saxena, K. Ahilan, J. Rivas, Phys. Rev. B 79, 060401 
(2009). 

M. Chattopadhyay, P. Arora, S. Roy, J. Phys.: Condens. Matter 21, 296003 (2009). 

N. Khan, A. Midya, K. Mydeen, P. Mandal, A. Loidl, D. Prabhakaran, Phys. Rev. B 82, 
064422 (2010). 

M. E. Fisher, S. K. Ma, B. G. Nickel, Phys. Rev. Lett. 29, 917 (1972). 

K. Ghosh, C. J. Lobb, R. L. Greene, S. G. Karabashev, D. A. Shulyatev, A. A. Arsenov, Y. 
Mukovskii, Phys. Rev. Lett. 81, 4740 (1998). 

S. F. Fischer, S. N. Kaul, H. Kronmuller, Phys. Rev. B 65, 064443 (2002). 

S. V. Grigoriev, S. V. Maleyev, A. I. Okorokov, Yu. O. Chetverikov, R. Georgii, P. Boni, D. 
Lamago, H. Eckerlebe, K. Pranzas, Phys. Rev. B 72, 134420 (2005). 


12 


201711.02794v1 


chinaXiv 


TABLE I: Comparison of critical exponents of MnSi with different theoretical models and related 


materials (MAP = modified Arrott plot; ND = neutron diffraction; Hall = Hall effect; AC = ac 


susceptibility; SC = single crystal; PC = polycrystal; cal = calculated). 


Composition technique Ref. To(K) B y ô 
Mnsi>® MAP This work 30.5  0.242+0.006 0.915+0.003 4.734+0.006 
Tricritical mean-field theory - 0.25 1.0 5.0 
Mean-field theory - 0.5 1.0 3.0 
3D-Heisenberg theory - 0.365 1.386 4.8 
3D-Ising theory = 0.325 1.24 4.82 
MnSi?? MAP 29.67 0.238+0.001 1.2040.01 6.140.1°" 
MnSi°C ND 28.82 0.22(1) - - 
Feo.sCo0.28i1P C Hall 36.0 0.371+0.001 1.380.002 4.78+0.01 
Cu20Se03° AC 58.3 0.37(1) 1.44(4) 4.9(1) 
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FIG. 1: (Color online) (a) The temperature dependence of magnetization [M(T)] for MnSi (the 
inset shows the magnification of the transition); (b) the isothermal magnetization [M(H)] at 4 K 


(the inset gives that in the lower field region). 
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FIG. 3: (Color online) The isotherms of M'ê vs (H/M)"/7 with (a) 3D-Heisenberg model; (b) 


3D-Ising model; (c) tricritical mean-field model; (d) the normalized slopes [NS = S(T)/S(Tc)] as 


a function of temperature. 
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FIG. 4: (Color online) (a) The temperature dependence of Mg and yo! for MnSi; (b) the isother- 


mal M(H) at To with the inset plane on the lg — lg scale (the solid curves are fitted). 
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FIG. 5: (Color online) Effective exponents (a) ßerp and (b) Yerr as a function of the reduced 


temperature £ for MnSi (dashed curves are guided on eye). 
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FIG. 6: (Color online) (a) Scaling plots of renormalized magnetization m vs renormalized field h 
around the critical temperatures for MnSi (only several typical curves are shown); (b) m? vs h/m 


below To in lower field region (the inset plot the Hc vs T with the linearly fitting solid line). 
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